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Abstract 



B° s — > £>,T and B® —> D*~ weak transition form factors are estimated for the whole physi- 
cal region with a method based on an instantaneous approximated Mandelstam formulation 
of transition matrix elements and the instantaneous Bethe-Salpeter equation. We apply the 
estimated form factors to branching ratios, CP asymmetries and polarization fractions of non- 
lcptonic decays within the factorization approximation. And we study the non-factorizablc 
effects and annihilation contributions with the perturbative QCD approach. The branch- 
ing ratios of semi-leptonic B® — » Di*^ l + ui decays are also evaluated. We show that the 
calculated decay rates agree well with the available experimental data. The longitudinal 
polarization fraction of B s — > D* S V{A) decays are ~ 0.8 when V(A) denotes a light meson, 
(■-<■) ' and are ~ 0.5 when V(A) denotes & D q (q = d, s) meson. 

o' 
o, 

■ 1 Introduction 

o 

In the past few years, charmless non-leptonic B s decays have been extensively studied [1], how- 
ever, the decays of B s to charmed particles are relatively less studied. Therefore, it is of urgent 
interest to put more attention on this topic. Semi-leptonic B® — > decays and non- 

leptonic B® — > Ds* X + two body decays (and their conjugated ones), where X + denotes a light 
meson or a D q (q = d, s) meson, can reveal useful information about the Cabibbo-Kobayashi- 
Maskawa (CKM) phases, the B° s - B® mixing parameters [2, 3], and the physics of CP vio- 
lations [4]. Studies on B s decays to charmed particles can be used to check the factorization 
hypothesis [5] and to search for physics beyond the Standard Model (SM) [6, 7, 8]. 

For weak decays, the transition form factors play an important role. Theoretically, to esti- 
mate the form factors of a relevant process, one has to rely on some non-perturbative approaches 

such as the Bethe-Salpeter (B-S) equation, quark models, QCD sum rules (QCDSR) and lattice 

(*) 

QCD. Turn to B s —> D s weak transitions, several works have been done: Early works, e.g. 
[9, 10], usually relied on the famous Bauer-Stech-Wirbel (BSW) model [11]. Authors of [12, 13] 
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adopted QCDSR for the calculation. In [14], the form factors are estimated within the covari- 
ant light-front quark model (CLFQM). Authors of [15] used the so called light cone sum rules 
(LCSR) to investigate form factors at large recoil, and heavy quark effective theory (HQET) 
to describe them at small recoil region. Each of these methods sketches one or another profile 
of non-perturbative QCD, and each has advantages as well as shortcomings. So it is worthy to 
estimate the B s — > form factors in another method which is based on the B-S equation [16] 
and the Mandelstam formulation [17] of the transition matrix element. To make predictions on 
non-leptonic decays, there is another task, that is how to evaluate the decay amplitudes with 
form factors available. It is well known that factorization approximation (FA) [11] has been 
extensively applied in non-leptonic weak decays and has been justified to success in explaining 
the branching ratios of several color-allowed B q decays [18]. The works mentioned before all 
adopted the FA to evaluate non-leptonic decay amplitudes. However, estimations based on the 
FA still suffer uncertainties from the non-factorizable effects and annihilation diagrams contri- 
butions, especially for the CP asymmetries (CPAs). Thus approaches beyond the FA are in 
need. Till now several approaches which can cover the non-factorizable effects have been devel- 
oped, such as the perturbative QCD (pQCD) approach [19], the QCD improved factorization 
(QCDF) approach [20] and SCET approach [21]. Studies on B s decays into charmed particles 
with the perturbative QCD approach have been carried out in [22]. In this work, we evaluate 
non-leptonic decay amplitudes under FA, as well as estimate non-factorizable and annihilation 
contributions in the pQCD approach. Besides these direct calculation in the FA or pQCD, the 
authors in [23] used 811(3)^ symmetry to estimate the widths of a class of two-body B s decays 
with the help of experimental data of corresponding B decays. There have been some studies 
on semi-leptonic B® — > l + ui decays with the approaches such as QCDSR [13], LCSR [15], 
CLFQM [14] and constituent quark meson model (CQM) [24]. 

On the experiment side, the world averaged branching fractions of some of the B® — > 
X + decay modes are already available [25]. Recently, Belle Collaboration has reported 
the observations of B® — > , D*~n + and p + decays and measurements of their 

branching fractions [26]. Measurements of branching ratio for B® — > Djl + ui+ "anything" are 
given to be (7.9 ± 2.4)% [25]. However, the exclusive semi-leptonic decay rates for B® — > 

(*)— 4- 

D. s l + vi processes have not been measured yet. It is expected that in near future more and 
more channels of B s decays will be precisely measured experimentally. 

In this paper, we estimate B® — > Di*^ form factors by calculating the corresponding tran- 
sition matrix elements in the instantaneous approximated Mandelstam formulation with the 
wave functions obtained from the Salpeter equation. The Salpeter equation are derived from 
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the B-S equation in instantaneous approximation [16]. The benefits of this method are its 
firm theoretical basis and its covering relativistic effects. It is well known that the B-S equa- 
tion is a relativistic two-body wave equation. With the form factors calculated, we study the 
semi-leptonic B® — > D^~l + vi decays and the two body non-leptonic B® — > X + decays. 

For non-leptonic decays, we estimate the branching ratios, CP asymmetries and polarization 
fractions of processes in the FA, and we also estimate the non-factorizable and annihilation 
contributions with the pQCD approach for B s — > PP(PV, VP, VV) processes. Comparing our 
results with different theoretical predictions and experiment data would enrich our knowledge 
of B s weak decays to charmed mesons. 

The paper is organized as follows: In Section 2, a brief review on the Salpeter equation 
and the instantaneous approximated Mandelstam formulation of B® — > transition matrix 

elements (from which the form factors are extracted) are presented. In Section 3, we illustrate 
how we evaluate the decay amplitudes with the form factors calculated and how to evaluate 
the decay rates, CPAs and polarization fractions. Section 4 is devoted to numerical results and 
discussions. 

2 Form Factors of — > transitions 

Form factors are the crucial elements of decay amplitudes. In order to estimate B® — > 

form factors, first we use the improved Salpeter method illustrated in [27, 28] to obtain the 

(*) 

wave functions of B s and D. s mesons. In these literatures the authors solved the full Salpter 
equations instead of only the positive energy part of the equation. Now we give a brief review 
on this method. Under instantaneous approximation, the well-known B-S equation 



can be deduced to be the full Salpeter equation, which equals to the following coupled equations 



Here x P (l) is the B-S wave function of the relevant bound state. P is the four momentum of the 
state and p\, p 2 , mi, m 2 are the momenta and constituent masses of the quark and anti-quark, 




(1) 



[16]: 



(m — uj\ — uj 2 )^ + (q p± ) = A f {q p± Hq p± ) a£ (q P± ) , 

(M + Wi +U 2 )if-~(q p± ) = -A^(q p± ) V (q p± )A 2 (q p± ) 
^-(q P± )=0, tp-+(q p± )=0. 



(2) 



respectively, q is the relative momentum q = a 2 pi — a\p 2 , where a\ = 



1)11+1712 



and a 2 = 



mi+m2 " 
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V(P,k,q) is the interaction kernel which can be written as V(k p ^,q p ^) under instantaneous 
approximation. In our notations, q p always denotes -7= and q Pt = q — ^-P. The definitions 



of 7/(g£ 



A!, and ^1(2) are gathered as follows: 



dq 



2vr 



p x P {q), 



Af = (-1+ Af = -L^o* T (-2+ £J], 

w 1 = Jml-ql, uj 2 



m 



(3) 



A± (2) satisfy A+ (2) + A 1(2) = A± 2) f A± (2) = A± (2) and A± 2) &A* 2) = 0. With these A ± , the 
wave function 99 can be decomposed into positive and negative projected wave functions 



^iH P J = Af (Q P± )^ p (Q P± )£*£(QrJ- 
In deriving the coupled equations (2), the decomposition of the Feynman propagator 



(4) 



Sjipj 



+ 



(5) 



(-iy +1 q p + ajM - ujj + ie (-iy +1 q p + ajM + ojj - ie 

where j = 1 for quark and j = 2 for anti-quark has been used. 

The wave functions relevant to B s — > transition have quantum numbers J p = 0~ (for 
B°,B® and Df) and l~(for D*^) and are written as [27, 29] 

P i P 



M 



<Po-(+)(q P± ) 
¥i-(-)(q P± ) = (q P± -£ 



— a 1 (q p J + a 2 (q p J + 1 ^a 3 (q P J+ ^ 



■a 4 (q f 



h( qp± ) + ^b 2 (q P± ) + ^h(q p± ) + ^h(q 

1 



75, (6) 
+ M fh(q p ) 



+ f Pb 6 (q p± ) + (4 P± f - q P± ■ e x )b 7 {q p J + — (P p 4 P± - Pq p± ■ e*)b 8 (q P J 

where ai(q P± ) and bi(q p ^) are wave functions of q p ; M is the mass of corresponding bound 
state; e A is the polarization vector for J p = V state. In numerical calculation, Cornell potential 
is chosen as the kernel, and the explicit formulation is (in the rest frame): 

V{q)=V s {q) + V v {q)l°®l0, 



2 a s (q) 



(7) 



3vr 2 (g 2 + a 2 ) ' 

where the QCD running coupling constant a s (q) = r,J^£ N , — , t^jti — r; the constants A, a, a, Vq 

00 zimj log(a+q /jVq Cd j 



and Aqcd are the parameters characterizing the potential, which are fixed by fitting the exper- 
imental mass spectra. The parameters used in this work are m& = 4.96 GeV, m c = 1.62 GeV, 
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m s = 0.5 GeV, A = 0.21 GeV 2 , a = 0.06 GeV, a = e = 2.7183, A QC d = 0.27 GeV and for 
0- state, V = -0.432 GeV (cs), -0.212 GeV (6s), for 1~ state V = -0.212 GeV. With these 
parameters, the wave functions of interested mesons can be obtained by solving the coupled 
equations (2). The details of how to solve the equation (2) could be found in [28]. 

We now turn to evaluate weak transition form factors. The starting point is the Mandelstam 
formulation of transition matrix elements. Since the relevant wave functions are obtained from 
the B-S equation with instantaneous kernel, instantaneous approximation should be applied to 
the Mandelstam formulation which has been carried out in details in [30, 31, 32]. We now follow 
[30] to sketch the derivation of the instantaneous Mandelstam formulation as follows. According 
to Mandelstam formalism, the transition matrix element between two bound states induced by 
a current T M , e.g. 7 M , 7^75, is written as 

(mm^m)) = j ^^Tr[xi t i qf )r^ p ^ (8) 

where q\ and q 2 are the relevant quark fields operators. Here and hereafter in this section the 
superscript or subscript i and / denote the quantities of the initial state and the final state, 
respectively, in the transition. Using the instantaneous B-S equation and decomposing the 
propagators into positive and negative parts (see equation (5)), equation (8) can be deduced to 

(f(P f )\(q 1 ^ q2 MP l ))= (9) 

v(qf Pf± )*t(qfp. ± )r M A|(% Pi± Hq ip . ± )A+ (g iPi± ) 



J (27T)4 ij 



(a 1 fPf Pi + q fp _ -Qi + ie)(a 1 M i + q ip . - w\ + ie)(a 2 Mi - q ip . -to 2 + it) 



+ 



where Af (g /p . ± ) = ^[ jJrWii^i+mi/)] andw 1(2) = ^ (2)J -p? (2) withp 1(2) = a 1{2)f P fPi± ± 
qf P± - The ". . . " represent the terms involving negative part (A~). Since the contributions of 
negative energy parts are small, we can ignored them in equation (9) [30]. Actually we com- 
pared tp~~ and M -Pf p T +uj 1 f+Lo 1 - ' wn i c h arise from the terms with A~('s), to tp ++ of B s meson 
numerically. The wave function of B s drops to ~ when \q\ ~ 2.1 GeV, so the part with 
\q\ > 2.1 GeV makes rare contributions. It is found that for the main part of the wave functions 

< 2.1 GeV), <p— and Mi _ P f p ^ lf+u)u are ~ - 2% of <p++ as \q) increases. 

After integrating equation (9) over qi p , one obtains 

(f(P f )\(q^q 2 ) m )) = J ^[^J^ + (^)|]. (10) 

where ^ = 7 o^ ++ 7o and L r = ^X^ Aj-( g/p<± ); q fPf± = q fp ^ - '-^j^Pf + 
s r(wi p i ~ lffi P f)i witn s r = a 2fPfP, ~ u 2 . In calculation, the relation 1 = jf-jf- = 

f J J 
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(Ai~(g/ P ) + Aj~ {<lfp f± ))jfj has been used and again the negative part is ignored. The in- 
stantaneous transition matrix element used here, i.e. equation (10), is different from the one 
in [31]. In [31], the instantaneous approximation is done in the initial particle's rest frame for 
both the initial particle and the final particle in the transition matrix element, whereas in this 
method, the instantaneous approximation is done in the relevant particle's own rest frame. It 
is found that, for the transitions of B c decaying to charmed particles, the two methods are 
consistent with each other; whereas for the transitions of B s decaying to charmed particles, the 
former generally gives larger form factors. 

Recent years another method, in which the B-S equation also plays an important role as in 
our method, was extensively studied [33] and applied to describe meson observables [34, 35]. 
In the method (referred as the DSE method for convenience), the rainbow-ladder truncation, 
which is a symmetry-preserving truncation that grantee the axial-vector vertices satisfying the 
Ward-Takahashi identity, is applied to the kernel of the Dyson-Schwinger equation for a quark, 
i.e. the gap equation, and the B-S equation's kernel. By solving the gap equation and the 
B-S equation for considered channels, interested observables could be estimated. The DSE 
method has impressed us due to the success on describing the pion as both a Goldstone mode, 
associated with dynamical chiral symmetry breaking (DCSB), and a bound state composed of 
constituent u- and d-quarks. The DSE method is mainly different from our method in two 
points. First, the quark propagators are different; and second the B-S kernel are different. The 
quark propagator in the DSE method, as the solution of the gap equation, is the dressed quark 
propagator which is characterized by a momentum-dependent mass function. In our method, 
we use free quark propagators with constituent quark masses. Of course, propagators appear 
in the B-S equation or considered amplitudes should be the dressed one, however, studies on 
the gap equation shown that the propagators of u, d, s quarks receive strong momentum- 
dependent corrections at infrared momenta while the mass functions in heavy quark b and c 
propagators can be approximated as a constant. Thus significant differences appear in the light 
quark propagators. The propagators in the DSE could also exhibit confinement characters of 
QCD. The B-S kernel in the DSE method and in our method both have the one gluon exchange 
interaction which is described by the products of the strong running coupling constant and the 
free gluon propagator. Except that we use an instantaneous kernel while the DSE method does 
not, the differences of the kernels are, our kernel also involves a confinement potential, while 
in the DSE method the confinement is described by using quark propagators with no Lehmann 
representation. The basic concepts of the DSE method is attractive, however for applications 
involving a widely ranges of observables, further assumptions and parameterization are usually 
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adopted [34]. For example, the B-S amplitude for a heavy meson used in [34] are obtained not 
by solving the B-S equation, but by assuming a parameterized form and then fitting the data 
to fix the parameters. Furthermore, the forms of the B-S amplitudes of heavy mesons are too 
simple compared to ours. Despite these differences in the propagators and the B-S amplitudes, 
the form of a transition matrix element in the impulse approximation in these works is the same 
as ours: the Mandelstam formulation. 

Due to the argument of Lorentz covariance, the transition matrix element can be decomposed 
into several parts, where the form factors show up. As usual, we denote form factors by the 
following decompositions: 

(DJ\W\B°) = f+(Q 2 )P» + f-(Q 2 )Q», (11) 
(D*-\V»\B° S ) ee -i M .l M MQ 2 )^' PtP f, (12) 

(D*-\A^\B° q ) ee A(Q 2 ) e *' Pi P M 

+h ^ M l+ P M f Q " + /o(Q 2 )(M, + M/)e* M , (13) 

where P = + P f and Q = P t - P f . f±(Q 2 ), fv{Q 2 ) and MQ 2 ) (i = 0, 1,2) are the form 
factors of weak transition B® — > D^*^ . 

3 Non-leptonic two-body decay rate and its CP asymmetry for 

>/;; ' x 

In this section, we first treat the non-leptonic two-body decay in the framework of the factoriza- 
tion approximation, and then introduce the pQCD method to estimate the contributions from 
the non-factorizable effects and the annihilation diagrams. For B® — > + L~ decays induced 

by b — > c transition, where L~ denotes a light meson, the low energy effective weak Hamiltonian 
is given by 

n cS = ^^{CiMQi + C 2 (fi)Q 2 ], (14) 

where Q\ = (c a b a )v-A{<lf3Up)v-A and Q 2 = {c a bp)v-A{Qpu a )v-A with q = d or s. And for the 
double charmed B® decays, the low energy effective weak Hamiltonian for the AB = 1 transition 
is [36], 

H cS (AB = l) = ^ Yl V^AC^Ql + CzitiQl + £Ci(M)Qi)+>i-c, (15) 

V 2 p=u,c { i=3 J 

where V pq is the CKM matrix element with (p = u,c) and (q = d,s). Cj(/x) are the Wilson 
coefficients. The local four-quark operators Qi can be categorized into three groups: the tree 
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operators in b — > p transition Q\, Q^, the QCD penguin operators Qi (i = 3,4,5,6), and 
the electroweak penguin operators Qi (i = 7,8,9, 10). All these local four-quark operators are 
written as 



(16) 



(17) 



Ql = (qaPa)v-A(ppbp) V -A, 

Q\ = {QocPp)v-A{ppb a )v- A , 

Q3 = (qaba)v-A^2(q x pQxl3)v-A, 
Qi = (q~ab/3)v-A'^2(qxl3qxa)v-A, 

Q5 = (qab a )v-A^2(q x pqxp)v+A, 

qx 

Qe = (q a bp)v-A^2(QxpQxa)v+A, 

q x 

Ql = 2 (qgba)v-A ^ e qx (q x pq xl 3)v + A, 
qx 

Qs = ni.qgbpjv-A e qx (q x pq xa )y + A, 
q x 

Q9 = -^(qabg)v-A e Qx {q x pqxp)v-A, 
q x 

QlO = -Aqab^V-A y eq x (qxl3qxa)v-A, 
q x 

where q x ranges from u, d, s to c. The subscripts a, j3 are color indices. The operator (^i a ip20)v±A 
= V'ia7 M (l ± 75)^2/3- As usual, we define the combinations ctj of Wilson coefficients 

£121-1 = ^21-1 + -^, 0,2i = C21 H (Ify 



(18) 



where N c is the number of quark colors and is taken as N c = 3. 

Unde 
as [8, 11] 



Under the FA, the matrix element of B® — > Df^X two body decays can be factorized 



A = (DW+\(cb) V -A\B o s )(X-\(qp) v _ A \0), (20) 

where (X~ (Px)\(qp)v-a\Q) = ifo±Px wnen the meson X denotes a scalar (pseudoscalar) , and 
(X~ (Px)\(qp)v-A\0) = ifi+Mxe*^ when A is a axial vector (vector). / ± and fi± are decay 
constants of particle A. The decay amplitudes for the B" — >■ D s +L~ decays can be expressed 



as 



M = ^=V cb V: qai A. (21) 



The decay amplitude of double charmed decay can be written as [8] 



M = ^{\ c a 1+ £ A P K + < + £« + <)]} A, (22) 



p=u,c 



where X p = V p bV* and a? = a,i + If with If given as follows: 



if = i? = ^{Ci[y-G(m p> fc 2 )]}, (23) 
if = ifo = g^K^ + C ^ c )[y " G(m p , k 2 )]]. (24) 

2^ 



The penguin loop integral function G(m p ,k ) is given by 



, f 1 m 2 . - fc 2 x(l - x) 
G(mp, fc ) = -4 / x(l - x)ln— ^ ^ -dx, (25) 



where the penguin momentum transfer k 2 = -^-(1 + (^L ~~ m q)(l ~ ~ 3 %~)/M x + ( m q + 2.,„ i; 
M x)/ m l) P 7 ] • The £ in equation (22) arises from the contribution of the right-handed currents 
and depends on the J p quantum numbers of the final state particles. The collected expressions 
of £ are shown as follows: 

im,t-. — m,„Mm„-\-m l ~') ' s V / 



{m b -m c )(m c +m q ) ' 
2M|. 

(m b ~m c )(m q -m c ) ' 
2M|. 



£> a X(0+) 

d;x(o-) (26) 



(m b +m c )(m c +m q ) ' 
+ ( mi ,+m c )K- m(; ) ' ^s X ( 0+ ) 

, D^l*) 

where X denotes aD ? (q = s,d) meson with its J p shown in the bracket just following it. The 
current quark masses encountered in G(m p ,k 2 ) and £ are taken from [25] and then evolved to 
the scale fx ~ mi by the renormalization group equation of the running quark masses [36] : 

mM= " ,w i^w/ i 1 + ( ?ft— — s — )' (27) 

where 

97 10 

7m0 = 6C F , 7mi = C F (3C F + -iV c - y /), (28) 

N 2 — l 

and CV = £ Nc . The number of quark flavors is denoted as /, which is taken as / = 5 in the 
present paper. 

As indicated before, under the FA non-factorizable effects and the annihilation contributions 
are both neglected. But for double charmed B s decays they may contribute conspicuously, espe- 
cially for CPAs, since it has been indicated that the annihilation diagrams usually make domain 
contribution on the strong phases according to the pQCD analysis. Thus in this work, we esti- 
mate the non-factorizable and annihilation contributions in the pQCD approach to make more 
reliable predictions on non-leptonic decays. Concretely, we will estimate the contributions from 
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non-factorizble and annihilation (if exits) for B® — > + X~ decays, where only pseudoscalar 

or vector present in the final state. Here we will not illustrate the pQCD approach in detail, 
instead we refer the readers to the original paper [19] and [22] which provides the specific pQCD 
studies on B® — > + X~ decays for details of this method. Thanks to the efforts did by the 

authors of [22], from which we borrow the expressions of decay amplitudes of non-factorizable 
and annihilation diagrams. But we use our form factors in estimating the factorizable color- 
favored diagrams' contributions. For example, the decay amplitude of the B® — > D~ decay 
could be expressed as 

M = ^{v cb V: d [F^(a 1 ) + F^(C l ) + S F ]+V ub V,: d S F )}, (29) 
S F = F e LL (a A + a l0 ) + F e L n L (C 3 + C 9 ) + F' SP ( a(i + as) + F e L n R (C 5 + C 7 ) 

+F a LL (a 4 - aio/2) + F^(C Z - C 9 /2) + F s a p (a 6 - a 8 /2) + F^ n R (C 5 - C 7 /2), 

where each F^ L a ^ P \- • • ) corresponds to a certain diagram's contribution. The subscript "e" 
represents factorizable emission (color- favored) diagrams; "en" represents non-factorizable emis- 
sion diagrams; "a" and "an" represents factorizable and non-factorizable annihilation diagrams 
respectively. The superscript "LL" , "LR" and "SP" correspond to the contributions from the 
(V-A)(V-A) operators, the (V-A)(V+A) operators and (S-P)(S+P) operators respectively. In 
this work, we calculate the factorizable emission contributions, which can be expressed in terms 
of form factors and decay constants, with our estimated form factors; and calculate the other 
contributions with the pQCD approach. The exact expressions of F e {n)^a(n) can be found in 
[22]. 

The decay width of a two-body decay is 

B s pol 

where p is the 3-momentum of one of the final state particles in the rest frame of B®. Another 
important physical observable is CP asymmetry. Generally the amplitude for double charmed B® 
decays considered here can be written as M = V c bV* q T\ + V u bV* q T2. The direct CP asymmetry 
arise from the interference between the two parts of the amplitude and is defined as 

T(B^f)-T(B^f) 



"" r(so->/) + r(BO->/) 

ej2sin<5sin7 



\Gi/G 2 \ + \G 2 /Gi\ + ei2cos5cos7 

= D 1 Sln7 , 
1 + u-i cos 7 



(31) 
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where the weak phase 7 = arg(— y*y d ) o± arg( y^y" 3 ), the strong phase <5 = arg(Ti) - arg(T2), 

cb °d cb cs 

d = V^H, G 2 = V r u6 V r u * ? T 2 and 

_ 6j2sin(5 _ 6j2cos(5 

1 = \Gi/G 2 \ + \G 2 /G l y 2= \G l /G 2 \ + \G 2 /G l y ( ' 

with e\ = +1 for q = s, and e 2 = — 1 for g = d, respectively. 

Besides the branching ratios and CP asymmetries, the polarization fraction of B s — > VV(A) 
decays is another important observable. To illustrate the polarization fraction, one can write 
the decay amplitude as [38] 

M = h±M v(A)2 [a P?P? + b<T + ic e u ^ P ^V vl ,{Xi)e v[A)2u {\ 2 ), (33) 

where Pi (Mi), Py\ (My±) and Py(A)2 iMv{A)2) are the momenta (masses) of the initial particle, 
the particle picking up the spectator quark in the final state and the other meson in the final state, 
respectively. Coefficients a, b and c are defined as a = 2C/i/(Mj + Myi), b = foC(Mi + Myi) 
and c = 2Cfy /(Mi + Myi), where C denotes the term involving coupling constant, relevant 
Wilson coefficients and CKM matrix elements in front of the hadron matrix element A. X's are 
the helicities of the final particles. Then the decay amplitude of various helicities can be given 
as 



aM?P vl + b(P vl + P vl P v{A)2 ) 



(34) 



Myi 

M\\=V2b h±M v(A)2 , M ± = V2c f\±My {A)2 Mi\P vl \, (35) 

where Ml, M\\ and M± denote longitudinal, transverse parallel and transverse perpendicular 
part of the amplitude, respectively. The expressions of .M's apply under the FA. For pQCD 
calculations certain terms corresponding to non-factorizable and annihilation diagrams should 
be added to each M. The momentum Py\ and energy Pyi(Py^ 2 ) are taken in the rest frame 
of the initial particle, i.e. B®. The polarization fraction is defined as Ri = \M L \ 2 +\M^\ i +\M±\' 2 ' 
where i = L, \\ and _L. 

4 Numerical results and discussions 

4.1 Form Factors of B° s — > Transition and Semi-leptonic decays 

By solving the Salpeter equation (2), we obtain the wave functions of B®, Df, D*/ 1 mesons. 
Then we calculate the form factors of B® — > transition in the whole physical region 

numerically with equation (10). In calculation, the particles' masses Mgo = 5366.3 MeV, 
M D ± = 1968.47 MeV and M D *± = 2112.3 MeV [25] are used. The results are drawn in Fig. 1 
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and Fig. 2 for D~ and D*~ , respectively. The parameter-dependent uncertainty can be estimated 
by varying the input parameters of our model in a reasonable range. In this work we vary the 
parameters m^, m c , m s , A and Aqcd by ±5% to give the errors. In this section only, we denote 
the momentum transfer q = Pi — Pf, instead of Q. It should be noticed that our theoretical 
estimations suffer other uncertainties arising from the instantaneous approximation, since the 
s-quark is not a heavy quark. It is well known that describing the inter-quark interactions by 
a QCD-inspired potential (which relates to the instantaneous kernel) works well for a meson 
consisting of a heavy quark and a heavy anti-quark (i.e. b(b), c(c)), but is questionable in 
describing light mesons. For the present case, as s-quark is not heavy enough, assuming the 
(anti-)quarks interact instantaneously may cause (maybe sizeable) uncertainties, which also 
means that retardation effects may give contributions (maybe sizeable). Till now the problem 
encountered here is not totally solved. But part of the retardation effects have been studied 
in [39]. In that work, the authors assume the confinement kernel as V s ~ l/(— <?o + an< ^ 

2 

approximate it to be -4(1 + Iq^/q 2 ) by expanding %. Then the go is replaced by its "on-shell" 
value, which are obtained by assuming that quarks are on their mass shells. The "on-shell" 
approximation imply that the considered meson should be a weak binding system. Finally, the 
total effect is adding a term w i^^. a iyi ( \j (q — k) 2 + m 2 — + m 2 ) 2 to V s . m and k are the 
constituent mass and momentum of the (anti-)quark. For heavy-light system, it is better to 
take m and k to be the heavy quark's mass and momentum. By using such a interaction kernel, 
some of the retardation effects could be incorporated in calculations. Of course, this method 
didn't solve the problem totally, because only some of the retardation effects are incorporated 
and we don't know how much they are. Thus we won't take this interaction kernel as a corrected 
version of our potential in equation (7). But in order to obtain a qualitative feeling about the 
uncertainties arising from instantaneous approximation, we use the interaction kernel presented 
in [39], which incorporated some retardation effects, to estimate the mass spectra and form 
factors and compare them to our results without retardation effects. It is found that the relative 

variations between the two sets of results are: A J^ Bb ~ 0.3%, A A f Ds< ~* ) ~ 5% and A /^% ~$ ~ 

m b s M Ds(*) im =0) 

4 — 7%. Due to the reason indicated before, we emphasize that the actual errors (caused by 
describing the inter-quark interaction with a potential) may be larger. 

In Table 1, we compare our form factors at q 2 = with those from other approaches. This 
can be seen from the table: for P>® — > D~ transition, our results are roughly consistent with 
the results of BSW model and QCDSR, but larger than those of LCSR, therefore, it is expected 
that the LCSR method may give smaller decay rates for Bg D~ + L + non-leptonic decays in 
which the momentum transfer is near q 2 = GeV 2 . For B® — > D*~ transition, our results are a 
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Figure 1: Form factors of B® — > D s weak transition. 
Table 1: Form factors of B° -»• Dj and £ s ° D*~ transitions at q 2 = GeV 2 . 





This work 


BSW [9] 


QCDSR [12] 


LCSR [15] 


CLFQM [14] 


/+(0) 


u -°'-0.03 


0.61 


0.7 ±0.1 


0.43 




/-(0) 


u - L '-0.04 






-0.17 




MO) 


u - ' u -0.04 


0.64 


0.63 ±0.05 




0.74 ±0.05 


/o(0) 


U - DO -0.01 


-0.56 


-0.62 ±0.01 




-0.61 ±0.03 


/i(0) 


u -°'-o.oi 


0.59 


0.75 ±0.07 




0.59 ± 0.04 


/ 2 (0) 


-0 70 +u - U3 











little larger than the results from other methods. 

For semi-leptonic B® — > l + ui decays induced by b — > cuj transition, the effective 

Hamiltonian can be written as [15] 

Gf . 



K 



eflf 



V2 



F c *6 7m (1- 75 )c^ 7 m(i_ 75 )/. 



(36) 



The amplitude of B® — > l + i>i decays could be obtained by sandwiching equation (36) 

between the initial and final states, which reads 



(37) 



The width of a semi-leptonic decay is T = f ]C po i \-M\ 2 dEidEf where E[ and Ej are 

the energies of the lepton and the meson in the final state respectively. With the form factors 
calculated, we estimate the branching ratios of semi-leptonic B® — > l + vi decays. The results 
are listed in Table 2 together with those from other approaches. The orders of magnitude for 
those branching ratios are consistent with each other. The B® — > e + v e decay rates have 

been studied in the same model in [40] as used here. But the parameters and the formulation 
of the transition matrix elements used in that reference are different from those in this work, 
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Figure 2: Form factors of B® — > D*~ weak transitions. 



Table 2: Branching ratios for semi-leptonic B® — > D s * l + i>i decays compared with those from 
other approaches. 



Processes 


This work 


CQM [24] 


QCDSR [13] 


LCSR [15] 


CLFQM [14] 


B»^D7l+v t {l = e,ix) 


(1.4 - 1.7) % 


(2.73 - 3.00)% 


(2.8 - 3.8)% 


i.oi-j % 




B° -> DjT+Vr 


(4.7- 5.5) x 1CT 3 






3.3+it x 10" a 




B°^-D* s -l+vi (Z = e, M ) 


(5.1 - 5.8) % 


(7.49 - 7.66)% 


(1.89 - 6.61)% 




5.2t£» % 


B° -> D*~t+u t 


(1.2 - 1.3) % 








i-si'o.i % 



thus their results are larger. The parameters used in calculation such as V c b, t b o are shown in 
the next subsection. 

4.2 Non-leptonic B s decay 

Now we can use the form factors to estimate the decay rates of B®. The CKM matrix elements 
used in our calculation are [25] 

\V ud \ = 0.97425, \V US \ = 0.2252, \V ub \ = 3.89 x 10~ 3 , 
|y cd | = 0.230, \V CS \ = 0.9735, \V cb \ = 0.0406. 
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Table 3: Decay constants used in our calculation in unit of MeV. 



u 


fk 


fD 


fD s 




fk* 


fa\ 


130 [25] 


156 [25] 


207^ [25] 


258^ [25] 


205^ [42] 


217± b 5 [42] 


229 [43] 


fo* 


I'd* 






/d s i(2460) 


/d s i(2536) 


f Bs 


245^ [44] 


272- ];; ii 


mtt [45] 


109l| [45] 


2271^ [45] 


77.3™ [45] 


0.23 [22] 



The lifetime t#o = 1.472 x 10 12 s [25] is taken in calculation. The Wilson coefficients are quoted 
from [41], and for fi ~ nib, they are: 

Ci = 1.0849, C 2 = -0.1902, C* 3 = 0.0148, C 4 = -0.0362, C 5 = 0.0088, 
C 6 = -0.0422, £t = -0.0007, = 0.0565, ^ = -1.3039, ^ = 0.2700, 

where a e (M\y) = 1/128 is the electromagnetic coupling constant. The strong coupling constant 
is taken as a s (mb) = 0.216. The decay constants used in this paper are shown in Table 3. 
Other inputs in pQCD analysis such as the wave functions and the Jet function appearing in 
the non-factorizable amplitudes and annihilation amplitudes are taken as the same as in [22]. 

With these input parameters, we calculate the branching ratios of non-leptonic B® to charmed 
particle decays. The results are listed in Table 4, together with those from other methods, as 
well as with available experimental data. The branching ratios of double charmed decays shown 
in the table are the CP averaged values: ^(B(Bg — > f)+B(B® — > /)), where B denotes branching 
ratio. The momentum transfer of B® — > L~ decays, where L~ is a light meson, is q 2 ~ 

— 1.5 GeV 2 , and the momentum transfer of the double charmed decays is q 2 ~ 3.5 — 6.4 GeV 2 . 

Compared our FA results and the results with corrections estimated by the pQCD approach 
(we will call them the pQCD corrected results), it can be found that the non-factorizable effects 
(and annihilation contributions when exist) for B® — > +L~ channels make little corrections 

to the FA results. For double charmed B s decays, the factorizable color-favored diagrams still 
dominate the decay widths, but contributions from other diagrams give up to ~ 20% corrections. 
We also estimate some B® — > + S(A) decay rates under the FA, where S represents scalar 
and A stands for axial vector, but the uncertainties due to non-factorizable effects should be 
realized. The results in the 2, 4-7 columns are calculated in the FA, but with different approaches 
to evaluate form factors. For B® — > D+X~ decays, the results from the FA are roughly consistent 
with each other except the values from light cone sum rules (LCSR) approach, which are smaller 
than others. This discrepancy reflects the difference of form factors, which has been shown 
explicitly in Table 1. For B® — > D* + X~ decays, our results are consistent with HQET and 
CLFQM values, but larger than those in BSW and QCDSR methods. 

Now we turn to compare the our results with the pQCD results in [22]. The major difference 
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Table 4: The (averaged) branching ratios of B® — > Ds* X~ in the unit of %. Our results 
with pQCD represent those with non-factorziable and annihilation contributions estimated in 
the pQCD approach. The parameters in [10, 12, 46] are replaced by the ones used in our paper: 
\V cb \ = 0.04, t b o = 1.472 x 10" 12 s, cti = 1.02. The results of B° s D+X' in the 7th column 
are quoted from [15] and the results of B® — > D* + X~ in the 7th column are quoted from [14]. 
Experimental data are quoted from [25, 26]. 



channels 


Ours (FA) 


Ours (pQCD) 


BSW |10] 


QCDSR |12] 


IIQET [46] 


[14] [15] 


pQCD [22] 


Experimental Data 


D + tt~ 


"• z '-0.03 


"■•^'-0.03 


0.37 


0.41 


0.31 


17+ u - uv 
ul '-0.06 


196+' llUli 


0.32 ± 0.05 


D t P~ 


64+ Uly 
U,D, -0.11 


60+ U11 
"'""-0.11 


0.94 


1.1 


0.82 


42+"- 17 

-0.14 


47 +0.24y 
"■^'-0.177 


± 0.11 ± 0.13 


D t k ~ 


"■ Uj!1 -0.002 


+ U.U01 

u.u^i_ 002 


0.028 


0.033 


0.023 


013+ UUUb 


0170+' UU)SV 

"■" ± '"-0.0066 




D t k *~ 


0.038_„ 005 


,, n , 8 +0.UU5 


0.049 


0.049 


0.040 




o.028i+|!;gi« 




°Ul 


n 7 c+u.ub 

-0.08 




0.88 


0.9 












o.3it«;S 


"• 3() -0.03 


0.28 


0.16 


0.31 


o.35±r i 


o.i89± ; 072 


0.24+ 1 '' u4 ' ± 0.03 ± 0.04 


D» 3 + p- 




°- 91+ o 


0.88 


1.1 


0.97 


iiS -0.31 




1.19+0 ± °' 17 ± °' 18 


D* + k~ 


"■"-^-0.002 


"■ uz ^_0.003 


0.020 


0.016 


0.023 


o n98+"- 00;i 

0.028_ 003 


0164+ IU)I)S4 

U.U1D«_ 0064 




DJ+fe* - 


°' 056 -0.007 


,, „ o + U.OUV 
U-Uo»_0 .007 


0.048 


0.049 


0.053 


,, ,,,,+0.006 
"'"""-0.006 


u-u*"_ .0124 




D a + a l 


1 3+ U1 




1.1 












D t D ~ 


4-lUJUB 
u.uoi_ 005 


"• 03(, -0.006 


0.027 


0.041 


0.049 


011 + "-'" 14 
"' uii -0.004 


022+' U)14 
"•""-o.oio 




D+D*- 


034 +0.009 


041+ U U11 
"'"* i -0.011 


0.031 


0.033 


0.034 


014 + ' UI0B 


021+ U Ula 
"'" zl -0.009 






012+° 

"'""-0.002 










002+ UUU1 
"'""^-0.001 






D* + Z}~ 


() 032+U.00B 


„ ,,.,,, + 0.006 

u.ud»_ 008 


0.013 


0.016 


0.036 


037+"- 004 

"'"•"-0.004 


027+ -" 17 

"•"'"-0.011 




D*+D*~ 


u - 11 -0.03 


13 +U.03 


0.074 


0.065 


0.11 


+ U.010 

0.086_ 000 


„ „+0.U26 
u '"' w -0.010 






o oio+ uuua 
















D t D 3 


n «3+ u lu 
u '*"-o.io 


1J -0.3 


0.51 


0.82 


1.2 


35+" 14 


,, + 0.36 
"'°°-0.24 


1.04 ± 0.35 




70+ U1(i 
"' '"-0.15 


1 0+ U ' 2 


0.57 


0.65 


0.77 


33+ U la 
u ' aa -0.11 


O dR+ U ' 31 
°' 48 -0.21 


2.75+g ± 0.69 


D t + D 7 


0.84_ 12 


1 2 +u ' J 

1 '^-0.1 


0.23 


0.33 


0.81 


92+"' 11 

"' az -0.11 


70+"- 44 
"■ '"-0.31 




2 4+ 11 - 4 

z '^-0.4 


JU -0.5 


1.48 


1.3 


3.2 


2 36+ 114 " 


99+ U r2 
" ,y -0.54 


., n „ + l.yy + 0.85 
o.uo_ 1 .04-0.86 


D f D "sO 


13+"" 1 
u ' lc> -0.02 










,, n _,+U.U22 
"'""' 1 -0.018 






0+0^(2460) 


37+ Uil! 
















D+ZJ- (2536) 


+u.uiy 

U.LMH_ ni3 


















i2+°- UJ 
ul ^-0.02 
















D* + D 3l (2460) 


1 ' 8 -0.4 
















£>* + 0^(2536) 


21 + "- ,l>< 
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between our pQCD corrected results and their fully pQCD results is that, for calculating factor- 
izable color-favored contributions, we estimate them in terms of our form factors, while in [22] 
all the contributions are calculated in the pQCD approach. Besides, some input parameters are 
different between the two works. The comparison can be seen from Table 4: most of the pQCD re- 
sults are smaller than our FA results as well as the pQCD corrected results, but close to LCSR re- 
sults. The Isgur-Wise form factor at maximum recoil in the pQCD is ^ s ~^ Ds = 0.44, whereas the 
corresponding quantity in our work is £ + = (Mb s (/++/_)+ Md 3 (/+ — /_)) /2 a/ Mb s Md s = 0.55. 
Besides, authors of [22] found that the form factors from their pQCD calculations are similar 
with those from LCSR. Therefore, we can simply state that form factors from the pQCD ap- 
proach in [22] are smaller than ours, which means that the discrepancies in decay rates between 
their pQCD results and ours could be attributed to the difference in form factors, rather than 
non-factorizable effects. 

Thanks to the efforts done by Belle, CDF, DO and other Collaborations, some experimental 
data of two-body non-leptonic B s decays are available now. We can see from Table 4 that our 
results agree well with the present experimental data. Since only some B" -»• DT' X~ decay 
modes have been observed and the experimental errors are still large till now, it is expected in 
the near future more precise tests could be made on theoretical predictions as the increase of 
the Bg events. 

(*) (*) 

We now investigate the direct CP asymmetries of B s — >■ Dl ' D\ ' decays. The results together 
with D\, L>2 are shown in Table 5. Di, D2 are defined in equation (32). Their values with the 
weak phase 7 determine the direct CP asymmetries. We show both the CPAs estimated under 
FA and in the pQCD approach. The discrepancy between FA and pQCD (or our pQCD corrected 
results) is obvious and is found mainly arising from the strong phases. It should be mentioned 
that another method called QCD improved factorization, which can cover the non-factorizable 
effects as well as pQCD, makes quite different predictions on CP asymmetries of some decay 
modes compared to pQCD. The reason is that the leading sources of the strong phase are 
different between the two approaches [47]. According to our results in Table 5, most of the 
direct CP asymmetries are too small to be tested experimentally for now. 

The results of polarization fractions of B® — > VV(A) decays are listed in Table 6, compared 
with other theoretical estimates and available experimental data. From the table, we can say 
that the non-factorizable contributions (and annihilation contributions when exit) do not give 
sufficient corrections on factorizations. Besides, it could be found that although our form factors 
are different from those in [22], the polariztions are still similar, which tells that the polariztions 
are less affceted by the form factors. Please note that for decays with a D* and a light meson 
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Table 5: The direct CP asymmetries in double charmed non-leptonic decays of B® in the unit of 
1CT 2 . Di and D2 are defined as Aq t p = D\ 1+ ^"cos7 ' wnere 7 is the weak phase (see equations 
(31,32)). In calculation, the weak phase is taken as 7 = 68.8°. We show the results with non- 
factorziable and annihilation contributions estimated in the pQCD approach. The FA results 
are also shown in the brackets. 



Final States 


Di in this work 


£>2 in this work 


Aq P in this work 


A^ P in [22] 


D+D- 


-0.36- -0.24 (5.3) 


9.8 (10) 


-0.3- -0.2 (4.8) 


n c+U. 0+0.1 
u -°-0.1-0.2 


DfD- 


0.12-0.15 (-0.32) 


-0.47 (-0.97) 


0.1 (-0.30) 




D+D*- 


-0.064- 0.039 (1.5) 


3.1 (2.6) 


-0.06- 0.04 (1.4) 


n -,+0.1+0.0 
u - 1 -o. 0-0.0 


DfD* s - 


0.048 - 0.063 (-0.081) 


-0.16 (-0.14) 


0.05 - 0.06 (-0.076) 




D*+D~ 


0.56 - 0.71 (-0.46) 


-0.69 (-1.6) 


0.5 - 0.7 (-0.43) 


u -*-o. 0-0.1 


D* s + D- 


-0.074- -0.064 (0.027) 


0.032 (0.13) 


-0.07- -0.06 (0.025) 




D*+D*~ 


0.-0.13 (1.5) 


2.9 (2.6) 


0. -0.1 (1.4) 


n i+O.O+O.O 
u - 1 -o. 1-0.1 


D*+D* s ~ 


-0.0018 - 0.015 (-0.081) 


-0.20 (-0.14) 


-0.002 - 0.01 (-0.076) 





in the final state, Rl ~ 0.8 S> R», R±; for decays with two charmed mesons in the final state, 
Rl ~ R\\ ~ 0.5 S> R±. The similar results have been found in B — > VV decays. 

Table 6: Polarization fractions of B® — > VV or V A decays. Rl, R\\ and Rj_ are longitudinal, 
transverse parallel and transverse perpendicular polarization fractions, respectively. Our results 
with pQCD represent those with non-factorziable and annihilation contributions estimated in 
the pQCD approach. 



Final States 




Ours (FA) 


Ours (pQCD) 


[22] 


Experiment [26] 


D*+ P - 


Rl 


U - 8 ' 4 -0.003 




0.87 


1.05i u ;f (stat)i u ;^(syst) 




R \\ 


U.1U4_ 004 


p n n+U.UUb 






D*+K*~ 


Rl 


U -841_o.005 


n QK7+U-UU3 
U.00 1 _o.003 


0.83 






R\\ 


n -1 00+0. 005 
u.ioo_ 006 


n n4 +0.004 
U.1U4_ 004 






D*+a^ 


Rl 


n 7QQ+0.UU5 
u. ( oo_ 007 










R\\ 


n 991 +U.UUY 
u - zzi -0.008 








D*+D*~ 


Rl 


U.OZO_Q QQg 


U.Q11_Q QQQ 


0.56 ±0.14 






R\\ 


U.410_Q QQg 


439+u.uuy 






D*+D* s - 


Rl 




n C07+0.006 


0.53 ±0.15 






R\\ 


0.439™ 


n 42 7 +0.006 
U - 4/ ' -0.012 






£> s *+£> s - 1 (2460) 


Rl 


436 +0 ' UU4 










R \\ 


u - oiz -0.008 








£> s *+£> s - 1 (2536) 


Rl 


U.4Z4_ 004 










Rn 


U.OZ 1 _ .007 









To conclude, we note that the form factors of B^ — > D~ and B® — > D*~ estimated in this 
work are roughly consistent with other theoretical estimates such as the QCD sum rule, the 
light-front quark model and so on. Using the derived form factors, branching ratios of semi- 
leptonic and non-leptonic B® decays to charmed particles are estimated. For non-leptonic decays 
we evaluate amplitudes under the factorization approximation as well as in the pQCD approach 
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(for non-factorizable and annihilation contributions). Our predictions on decay rates are found to 
agree well with the available experimental data. Direct CP asymmetries estimated in FA and in 
pQCD are quite different (even have different sign). Generally, the CPAs of B® — > D q are 
less than 1%. Polarization fractions of B s — > D*V(A) decays follow the similar rule as B — > VV 
decays, that is, for decays with a D* and a light meson in the final state, Rl ~ 0.8 3> R\\, R±, 
for decays with two charmed mesons in the final state, Rl ~ R\\ ~ 0.5 3> R±. 
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